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ABSTRACT 

Measurements of the frequency with which short-period planets occur around main sequence stars allows a direct prediction of the 
number and types of such planets that will be amenable to characterization by high-contrast instruments on future giant segmented- 
mirror telescopes (GSMTs). Adopting conservative assumptions, I predict of order 10 planets with radii R P = 1—8 and equilibrium 
temperatures < 400 K should be accessible around stars within 8 pc of the Sun. These numbers are roughly the same for both near- 
infrared observations of scattered starlight and mid-infrared observations of planetary thermal emission, with the latter observations 
demonstrating greater relative sensitivity to smaller and cooler planets. Adopting the conservative assumption that planets with R P = 
1 — 2 Rg, and 2—4 R® occur with equal frequency, I predict a 40% chance that a planet with R P = 1—2 Rg, and equilibrium temperature 
200-250 K will accessible to high-contrast thermal infrared characterization; this would be a compelling object for further study. To 
validate these predictions, more detailed analyses are needed of the occurrence frequencies of low-mass planets around M dwarfs, both 
in the Kepler field and in the solar neighborhood. Several planets already discovered by radial velocity surveys will be accessible to 
near-infrared high-contrast GSMT observations, including the low-mass planets a Cen Bb and (depending on their albedos) GJ 139c 
and d, GJ 876b and c, and r Cet b, c, and d; r Cet f would be amenable to thermal infrared characterization. Further efforts to model 
the near-infrared reflectance and mid-infrared emission of these and other short-period planets are clearly warranted, and will pave 
the way for the interpretation of future high-contrast characterization of a variety of planets around the nearest stars. 

Key words. Planets and satellites: general — Planets and satellites: detection — Instrumentation: adaptive optics — Techniques: 
high angular resolution — Methods: numerical 



1. Introduction 

A 1.1 Earth-mass (M©) planet candidate has been announced 
around a Centauri B (Dumusque et al., 2012), and a family 
of five candidates with minimum masses 2-10 M© has been 
announced around r Ceti (Tuomi et al., 2012). If confirmed, 
the discovery of so many low-mass planets around stars very 
near the Sun would be a strong sign that such planets are com- 
mon. Such a conclusion would be consistent with early results 
from the Kepler mission, which indicate that the frequency of 
planet candidates increases sharply toward smaller planetary 
radii (R P < AR®) for orbital periods P < 250 d (Howard et al., 
2012; Dong & Zhu, 2012; Fressin et al., 2013). These Kepler re- 
sults also suggest a dramatic increase in the frequency of small 
planets (radii Rp = 2—4/?©) orbiting cooler stars (Howard et al., 
2012). Because such stars are more numerous than sun-like stars, 
there should be many such planets in the Solar neighborhood. 
Because most of these planets orbit M dwarfs, their equilibrium 
temperatures T eq are quite cool: typically 200^1-00 K. 

The prospect of numerous temperate, small, low-mass plan- 
ets orbiting nearby stars is an exciting one. Such planets may 
lie in their host stars' "habitable zones" and be capable of sup- 
porting liquid water on any solid surface (Huang, 1959; Kasting 
et al., 1993; Pierrehumbert & Gaidos, 2011). These planets 
around nearby stars are also the only ones for which in situ 
exploration could be even remotely feasible in the foreseeable 
future (e.g., Crawford, 2011, and references therein). For now, 
the proximity of nearby systems and their attendant high photon 



fluxes make these planets the most amenable to remote charac- 
terization of their atmospheric and/or surface compositions. 

In recent years much effort has been devoted to preparing 
for observations of potential Earth analogues transiting cool M 
dwarf using, e.g., JWST (Deming et al., 2009; Kaltenegger & 
Traub, 2009). However, because even short period planets are 
unlikely to transit, the nearest accessible transiting planets will 
be around rather fainter, more distant stars. The durations of 
these transits will be short and their signal to noise low, so 
tens or hundreds of observations will be necessary for signifi- 
cant atmospheric measurements. Phase curve observations pro- 
vide higher duty cycles and may be able to characterize small, 
non-transiting planets around stars closer to the Sun (Selsis et al., 
2011; Maurin et al., 2012). Spacecraft capable of making such 
measurements have been proposed (e.g., Vasisht et al., 2008; 
Swain, 2010; Tinetti et al., 2012), but these observations will 
require extremely stringent spectrophotometric stability (better 
than 10 -4 over periods of days). Considering these various limi- 
tations, it is reasonable to explore alternative modes of exoplan- 
etary characterization. 

In this paper I present the first quantitative estimates of 
the local exoplanet population accessible to ground-based high- 
contrast observations, and a discussion of how such observations 
can directly constrain the demographics of the local planet pop- 
ulation. Though it has been appreciated for at least a decade 
that the next generation of ground-based giant segmented- 
mirror telescopes (GSMTs) can characterize small, cool plan- 
ets via such observations (Angel, 2003; Hawarden et al., 2003), 
Kepler's ability to measure planet frequency as a function of P, 
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Rp, and stellar effective temperature r e ff with unprecedented pre- 
cision permits the first informed estimate of the number of short- 
period planets that will be accessible to future large-aperture fa- 
cilities. 

Sec. 2 discusses simulations of statistically representative 
planetary systems around nearby stars and an assessment of 
high-contrast GSMT instruments' ability to characterize these 
planets. Sec. 3 then presents the number of planets accessible to 
such characterization under various assumptions of planetary oc- 
currence frequency and instrumental capabilities, and discusses 
the consequences of these assumptions. The main results are 
summarized in Table 1 . In brief, of order 10 planets are predicted 
to be amenable to such characterization at either near-infrared 
(NIR, < 2.5 jum) or mid-infrared (MIR, ~ 10/i»i) wavelengths, 
and this number varies by no more than a factor of 2 under rather 
more optimistic or pessimistic assumptions. In addition, there is 
significant likelihood (at least 40%, and probably ~4x higher) 
that such observations will be able to characterize an Earth or 
Venus analogue {Rp = 1 — 2R m , T eq = 200 - 250 K). Finally, 
I summarize the results of this work and suggest useful future 
work in Sec. 4. 



2. Methods 

2.1. Defining the Local, Observable Stellar Population 

To construct an observable sample of nearby stars, I take the 
8pc sample of Kirkpatrick et al. (2012) and remove all white 
dwarfs and spectral types later than M7 (to avoid brown dwarfs). 
White and brown dwarfs may have low mass planets, and these 
planets could even be temperate (Andreeshchev & Scalo, 2004; 
Monteiro, 2010), but at present the frequency of planets in such 
systems remains very poorly constrained. Because high-contrast 
observations have difficulty attaining the contrast necessary for 
planet characterization in the presence of nearby bright point 
sources (mitigation strategies exist but require customized in- 
struments; see e.g., Crepp et al., 2010; Cady et al., 2011), I re- 
move all known binaries with apparent separations < 5". Stars 
with fainter companions (e.g., Procyon A and Sirius A) are re- 
tained. For objects listed by Kirkpatrick et al. (2012) with a range 
of spectral types, I take the mean of the specified range. Some M 
dwarfs are listed without any spectral type, and for these I take 
the spectral types specified by Rojas-Ayala et al. (2012). The fi- 
nal list contains 92, 20, 6, 2, and 4 stars of types M, K, G, F, and 
A, respectively. 

Accurate stellar parameters are not known for all these ob- 
jects, so I adopt a homogeneous approach and assign T e g and 
stellar radii using empirical, interferometrically-determined val- 
ues (Boyajian et al., 2012a,b) for each spectral type. Apparent 
V-band magnitudes are taken from the UCAC4 Catalogue 
(Zacharias et al., 2012), converted to absolute V magnitudes us- 
ing the parallax values in Kirkpatrick et al. (2012), and con- 
verted to stellar masses using the mass-luminosity relationship 
of Henry & McCarthy (1993); using the K band relationship of 
Delfosse et al. (2000) affects the masses thus derived by roughly 
15%, but this does not significantly affect my results. Finally, a 
BT-Settl stellar model (Allard et al., 201 1) is generated by inter- 
polating in logg and r e g to the stellar parameters, and account- 
ing for the star's size and distance from the Sun. 

Note that this 8 pc sample is sufficiently small that high- 
quality spectra could be obtained for all these targets with rela- 
tive ease. Spectra at low resolution would provide a useful abso- 
lute flux calibration, and spectra at high resolution would permit 
stellar parameters and abundances to be derived at high preci- 



sion. Such a homogeneous stellar characterization effort would 
be highly desirable, and would reduce the need for the assump- 
tions about stellar parameters described above. 

2.2. Modeling the Local Exoplanet Population 

I estimate the distribution of the local exoplanet population us- 
ing an analysis of planet frequency as a function of r e ff, P, and 
Rp from the Kepler planet candidate sample derived from the 
first three quarters of Kepler observations (Bomcki et al., 201 1; 
Howard et al., 2012). Although these frequencies are computed 
for planet candidates and not fully validated planets, the occur- 
rence of false positives in the sample is below 15% (Morton & 
Johnson, 2011; Fressin et al., 2013). Furthermore, planet candi- 
dates with radii Rp = 2 — 4R<$ (the targets most accessible to 
high-contrast characterization) have a contamination fraction of 
< 7% (Fressin et al., 2013). Thus the Kepler false positive rate 
has at most a minor effect on the results presented here. 

Hot Jupiters are only half as frequent in the Kepler sample 
as observed by radial velocity surveys (Howard et al., 2012). 
The authors attribute this discrepancy to the lower metallicity of 
the Kepler sample, under which conditions fewer massive plan- 
ets should be expected (Santos et al., 2004; Fischer & Valenti, 
2005). If this discrepancy is confirmed, then the simulations pre- 
sented here underestimate the frequency of large planets by as 
much as a factor of 2; these planets are intrinsically rare, so this 
is a small effect. The frequency of Neptune-mass and smaller 
planets depends less strongly on stellar metallicity than does the 
frequency of Jupiter-mass objects (Sousa et al., 2008; Buchhave 
et al., 2012); since I show below that smaller planets comprise 
the bulk of the planet population accessible to high-contrast 
characterization, differences between the Kepler and local stellar 
populations should not strongly affect these results (this assump- 
tion will be tested in the future by high-precision radial velocity 
surveys). 

I begin by determining the frequency of planets of given Rp 
for each star, based on its r ef f. I assume that binary stars have 
planet frequencies consistent with those of single stars, and that 
4 — 8 /?© and 8 — 32 /?© planets occur with frequencies of 1 .7% 
and 0.79%, respectively (see Figure 8 of Howard et al., 2012). 
For 2 — 4/?® planets, I compute a frequency using Eq. 9 of 
Howard et al. (2012), which was defined using T e g in the range 
3600-7100 K. I extend this relation to lower and higher T e ff by 
assuming that beyond the specified range the frequency distribu- 
tion of planets is uniform and everywhere continuous. In Sec. 3.2 
I quantify the effect of this assumption. Howard et al. (2012)'s 
planet statistics are incomplete for Rp < 2/?®, but their results 
suggest 1 — 2/?© planets are at least as common as 2 — 4R® 
planets. I conservatively assume that planets in these two size 
ranges occur with equal frequency, consistent with more recent 
analyses (Dong & Zhu, 2012; Fressin et al., 2013). Throughout 
this work, I assume a uniform distribution of planet sizes within 
each of the Rp ranges noted above. 

The preceding steps determine the number of planets of vary- 
ing Rp expected around a given star for P < 50 d. In each of 
10 4 trials I assign to each star either no planet or a planet in 
one of the four Rp bins listed above, with the appropriate fre- 
quencies. I then assume that P in each radius range follows that 
given by Eq. 8 and Table 5 of Howard et al. (2012). Though their 
planet population statistics are limited to P < 50 d, even over 
this small range the frequency of planets increases in In P space 
for P > 10 d. I conservatively assume that planet frequency is 
flat for P > 50 d (see also Dong & Zhu, 2012). Letting planet 
frequency increase toward longer periods only slightly increases 
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the number of accessible planets, because the planets most ac- 
cessible to high-contrast characterization lie on short-period or- 
bits (see Sec. 3.2). Future studies should characterize planet fre- 
quency at longer orbital periods (this is a primary goal of the 
Kepler mission) and should explore the dependence of planet 
frequency on T e ff in greater detail. 

The result of this process is a planet population that pro- 
vides a direct estimate of the frequency of planets with Rp and 
P around each star in the 8 pc sample. Planetary systems may 
be more diverse than predicted by this simple model: for ex- 
ample, the M dwarf GJ 876 has multiple massive planets (e.g., 
Baluev, 201 1) while short-period planets of even very low mass 
have been excluded around the low-mass Barnard's Star and 
Proxima Centauri (Endl & Kiirster, 2008; Zechmeister et al., 
2009; Choi et al., 2012). Nonetheless the modeling approach 
described above has the advantages of being relatively simple 
and of providing estimates across a broad range of planetary and 
stellar parameter space that decently approximate the planet de- 
mographics presented by Howard et al. (2012). Future efforts 
should strive to quantify planet frequency over as broad a range 
as possible of T e ff, P, and Rp to obviate the need for some of the 
assumptions made above. 

2.3. Estimating High-Contrast Instrument Performance 

Adaptive optics (AO) systems designed to provide very high im- 
age quality are often termed "high contrast" systems, a name 
which emphasizes their sensitivity to faint objects located near 
much brighter ones; see Oppenheimer & Hinkley (2009) and 
Mawet et al. (2012) for recent reviews on this topic. At least one 
high-contrast system is already operating on-sky at the 5 m Hale 
telescope and reaching contrast levels of ~ 10~ 7 (Oppenheimer 
et al., 2012), and even more capable instruments are under con- 
struction for existing 8 m telescopes (Macintosh et al., 2006a; 
Beuzit et al., 2006). End-to-end simulations and testing of these 
latter instruments suggest they may reach final contrasts of 
roughly 10" 8 (Mesa et al., 2011; Thomas et al., 2011; Wildi 
et al., 201 1). At this level of performance, short-period exoplan- 
ets will likely remain beyond the reach of such instruments (see 
Sec. 3.3.1). 

The next generation of GSMTs will have diameters rang- 
ing from 25 m (GMT; Johns, 2008) to 30 m (TMT; Nelson & 
Sanders, 2008) to 39 m (E-ELT; Gilmozzi & Spyromilio, 2007), 
and all are projected to rely heavily on AO-assisted, diffraction- 
limited observations to maximize the science gains of these new 
facilities. High-contrast NIR instruments have been proposed 
and studied for all these telescopes (Macintosh et al., 2006b; 
Matsuo & Tamura, 2010; Brandl et al., 2010; Kasper et al., 
2010). These initial studies suggest that with moderate and rea- 
sonable technological advances, such instruments will be sensi- 
tive to objects at separations and contrast levels unprecedented 
for optical/infrared telescopes. 

High contrast instruments generally require highly special- 
ized designs to achieve their goals, and their performance varies 
across different targets, telescopes, and wavelengths. I estimate 
the performance of a generic GSMT high-contrast instrument 
using the analytic relations of Guyon (2005) and adopting the 
parameters in Table 4 of that work unless noted otherwise (vari- 
ations of instrument parameters are discussed in Sec. 3.3) 1 . The 
instrument is assumed to correct both phase and amplitude er- 



1 Note that Eqs. 16 and 17 of Guyon (2005) contain an error: the 
coefficient 0.484 should be 0.0484 (O. Guyon, private communication, 
Jan 2013). 



rors, sensed via a Mach-Zehnder pupil plane wavefront sensor. 
The wavefront sensing and science wavelengths are assumed 
to be identical to minimize chromatic contrast errors (Guyon, 
2005). The relations used compute raw contrast and not the ef- 
fective contrast limits for detection; I therefore assume that ob- 
servational and post-processing techniques can suppress the in- 
stantaneous PSF contrast by a factor of 30 at all angular sepa- 
rations; current techniques achieve suppression factors of 20 to 
several hundred at separations as small as 4.5A/D (Marois et al., 
2008; Crepp et al., 201 1; Vogt et al., 201 1), so the value assumed 
here may be somewhat conservative. Finally, following Lunine 
et al. (2008) I impose an absolute contrast floor of 5 x 10 -9 
in all cases. The result is a simulation of high-contrast perfor- 
mance with the phenomenologically correct dependencies on 
wavelength, telescope size, stellar flux, and other parameters 2 . 

In these performance estimates I also crudely account for 
the effect of finite stellar sizes, which are significant for sev- 
eral of the nearest and largest stars. For example, both a Cen A 
and Sirius A have angular diameters 0=8 mas, comparable to 
the diffraction-limited angular resolution (A/D) of a 30 m tele- 
scope observing at 1.2/im. Resolved sources cause the effective 
throughput of high-contrast observations to decrease markedly 
at the smallest separations (Guyon et al., 2006), though the 
attainable contrast is not adversely affected for proper choice 
of coronagraph or other starlight suppression system (Crepp 
et al., 2009). I approximate the deleterious effects of finite stel- 
lar size by imposing a minimum inner working angle (IWA) of 
(2.2 + 0.4 log l0 [6(A/ D)~ l ])A/ D, a relation which approximates 
the 50% throughput limit of an ideal coronagraph operating at 
10 -10 contrast (Guyon et al., 2006). This minimum separation 
increases for larger 8, more extreme contrast levels, and less op- 
timal coronagraph designs. In this work finite stellar size only 
affects a few star systems (those with the largest ff), and only in 
the 2A/D IWA case described in Sec. 3.3.1. Future high-contrast 
observations probing near or within IA/D (Serabyn et al., 2010) 
should consider these effects in greater detail. 

2.4. Estimating the Observability of Short-Period Planets 

I assess the accessibility of each simulated planet to high- 
contrast observations by considering both scattered starlight and 
thermal emission from the planet. I assume that all planets will 
be observed at quadrature, since radial velocity measurements 
will allow optimal phasing of most interesting targets. I fur- 
ther assume that all planets are on circular orbits and exhibit 
Lambertian scattering profiles, and that the observatory can ob- 
serve all stars in the target sample (under currently published 
plans GSMTs will be located in both hemispheres). 

I assign a Bond and wavelength-independent geometric 
albedo (A B and A ? , respectively) to each planet by drawing a 
value from a uniform random distribution spanning 0.0-0.4, a 
range which is a roughly representative for solar system planets 
(see also Hu et al., 2012). Although A B and A g describe different 
physical quantities (A g determines the strength of scattering at 
a given wavelength, while Ag determines a planet's equilibrium 
temperature r eq , and thus its thermal blackbody emission; see 
e.g., Seager, 2010), for simplicity I use the same value for each 
planet's A B and A g . This equality holds to within a factor of two 
for planets in the Solar System. At NIR wavelengths most ob- 
servable planets in the simulations below are only accessible in 



2 A Python software module used to calculate these quantities is 
available from the author's website, located (as of publication) at http : 
/ /www . mpi a . de /home s/ i anc / 
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reflected light, so higher or lower albedos will respectively in- 
crease or decrease the size of the accessible planet population; 
the effect is reversed for observations at thermal infrared wave- 
lengths. To each planet I also assign a heat redistribution factor 
(/; Seager, 2010) drawn from a uniform random distribution 
spanning the range (0.25, 0.66), where the extrema of this range 
respectively indicate zero and full redistribution of the incident 
stellar flux. The redistribution factor does not significantly affect 
the detectability of planets in reflected light, but would affect the 
thermal infrared observations described in Sec. 3.3.2. 

Using Ab, /, and the physical system parameters I then com- 
pute T eq for each planet, and I assume that the planet radiates as 
a blackbody with this temperature. Finally, the planet to star flux 
contrast in a given bandpass is the ratio of the sum of the planet's 
reflected and thermal photon flux to the photon flux from the 
host star. The planet is assumed to be observable if its contrast 
lies above the sensitivity limits computed in Sec. 2.3. 

3. Simulation Results 

In this section I discuss the local population of planets accessible 
to characterization by high-contrast observations under a num- 
ber of different assumptions about planet population parameters 
and instrumental configurations. I begin by presenting a set of 
baseline astrophysical and observational parameters in Sec. 3.1. 
Sec. 3.2 discusses the effects of varying the parameters of the 
local planetary population, and Sec. 3.3 discusses the impact of 
various telescope and instrument design choices. For all cases 
discussed, the number of planets of various Rp and T eq accessi- 
ble to high-contrast characterization is listed in Table 1 . 



3.1. Baseline Scenario 

The baseline observational parameters are a 30 m telescope ca- 
pable of observing down to an IWA of 3A/D: this instrument is 
analogous to TMT's PFI or SEIT instruments (Macintosh et al., 
2006b; Matsuo & Tamura, 2010), but with better wavefront sen- 
sor performance. This baseline case provides an angular resolu- 
tion midway between the 25 m GMT and 39 m E-ELT. I adopt a 
baseline observing wavelength of 1.2 fim (J band), where future 
AO systems should achieve good correction and which provides 
good angular resolution. Observations at J band are also inter- 
esting because of their potential to descry the 1.27 yum O2 band 
expected in oxygen-rich planetary atmospheres (e.g., Turnbull 
et al., 2006; Kawahara et al., 2012). In terms of astrophysical 
parameters, the baseline scenario assumes that 1 — 2/?© and 
2 — 4/?© planets occur with equal frequencies (as described 
in Sec. 2.2, this may be a conservative underestimate), and that 
planet frequency is flat in In P space for P > 50 d (the longest 
periods considered by Howard et al., 2012). 

Fig. 1 shows a representative planet sample and the com- 
puted sensitivity limits as applied to a Cen B. It is notable that 
the recently discovered low-mass planet candidate a Cen Bb 
(Dumusque et al., 2012) presents a contrast of 10~ 7 and is emi- 
nently observable in this baseline scenario (assuming a radius of 
1 . 1/?©, obtained by converting m sin i to R P via Eq. 1 of Lissauer 
et al., 2011). A simulation of the GJ 876 system is shown in 
Fig. 2; despite the worse contrast achieved on this fainter star, 
the massive planets GJ 876b and c (Delfosse et al., 1998; Marcy 
et al., 1998, 2001; Baluev, 2011) can be directly characterized 
if A g (Rp/R®) 2 > 19 and 6, respectively (A g > 0.02, with radii 
estimated as described previously). Though not shown, the cur- 
rently known planets GJ 139c and d (Pepe et al., 201 1) can also 
be characterized if A g (Rp/R®) 2 > 0.9 and 1.8, respectively (i.e. 
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Fig. 1 Simulated planet populations for a Cen B in the baseline 
scenario outlined in Sec. 3.1, observing at 1.2/j.m with an in- 
strument similar to TMT/PFI or TMT/SEIT. Planets above and 
to the right of the dashed line would be accessible to high- 
contrast spectroscopic characterization. The star symbol indi- 
cates a Cen Bb (assuming a radius of 1.1/?©). Each circle's size 
and color refers to a simulated planet's R P and r eq , as indicated 
in the legend. Most of the simulated planets shown are detected 
in scattered starlight, though thermal emission contributes to the 
shortest-period planets (including a Cen Bb). 
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Fig. 2 Same as Fig. 1, but for the M dwarf GJ 876. From left to 
right, the star symbols and error bars show the contrast expected 
for planets c, b, and e over the range of albedos and recirculation 
parameters described in Sec. 2.4. Planets b and c will be accessi- 
ble to characterization in scattered light if A g (Rp/R<s) 2 > 19 and 
6, respectively. 



A g > 0.4), and the planet candidates t Cet b, c, and d could be 
characterized for A g (R P /R®) 2 > 0.08, 0.4, and 1.2, respectively 
(A g > 0.04, 0.14, and 0.35). Theoretical efforts to model all these 
planets' NIR reflectance are clearly warranted. 

The number of accessible planets in this baseline scenario 
are shown in Fig. 3 as a function of T eq and Rp, and these re- 
sults are summarized in the first row of Table 1 . In all results I 
sort the planets by size (in the ranges described previously) and 
according to whether T eq is below, within, or above the range 
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Fig. 3 Predicted number of planets (as a function of Rp and r eq ) 
expected to be characterized in the baseline high-contrast GSMT 
observing scenario outlined in Sec. 3.1. The colors and the small 
numbers in each cell both indicate the expectation value for 
that combination of Rp and jT eq , and the histograms show the 
marginalized distributions. Note that the baseline scenario simu- 
lated here conservatively assumes that planets in the two smallest 
Rp ranges occur with equal frequency (see Sees. 2.2 and 3.2). 



appears that this relation does not apply to Rp < 2R®; nonethe- 
less, Table 1 shows that approximately equal numbers of plan- 
ets in the two smallest size ranges would be accessible to char- 
acterization if the trend held to smaller radii (i.e., if Kepler's 
completeness is significantly overestimated for Rp < 2/?©; 
Howard et al., 2012; Dong & Zhu, 2012; Fressin et al., 2013). 
However, because the smaller planets must be closer to their 
host stars to present the same contrast ratio as the larger planets, 
the additional 1—2 /?© planets have rather warmer temperatures 
(T^eq ~ 400 K) than the baseline population. 

My baseline scenario conservatively assumes that planet fre- 
quency flattens out for P > 50 d, consistent with recent results 
using additional Kepler data (Dong & Zhu, 2012). A recent anal- 
ysis of the first six quarters of Kepler data extends planet fre- 
quency estimates (Dong & Zhu, 2012). That work indicates that 
the frequency of planets in In P space for P < 250 d appears 
fiat for Rp < 4/?©, but continues to increase for larger plan- 
ets; this is consistent with radial velocity surveys, which show 
that the frequency of massive planets continues to increase out 
to 2000 d (Cumming et al., 2008). I investigate the impact more 
frequent planet occurrence on longer periods by letting Howard 
et al. (2012)'s planet frequency relations increase at longer peri- 
ods (instead of arbitrarily flattening the distributions). Table 1 
shows that compared to the baseline scenario, two additional 
planets (one each with sizes 4 — 8/?© and 8 — 32/?©) would 
be accessible to high-contrast characterization in this scenario. 
Because these additional planets occur on longer periods, they 
are rather cooler (r eq < 200 K) than the baseline population. 



200-300 K (a range which corresponds crudely to the T eq of 
Mars and a low-albedo Venus). Fig. 3 shows that the planets 
detected in the baseline scenario typically have T eq < 400 K 
and Rp < 4/?©: these planets are visible mainly via scattered 
starlight and not intrinsic thermal emission. Although Howard 
et al. (2012) show that planets with Rp < 4/?© are most fre- 
quently found around cooler stars, I find that few planets or- 
biting M dwarfs are accessible at distances >3.5 pc as a result 
of these stars' intrinsic faintness, and the consequent degrada- 
tion of achievable contrast. Finally, I estimate a 6% chance that 
in this baseline scenario an Earth/Venus analogue (1 — 2R®, 
T eq = 200 — 250 K) can be observed. This probability is likely 
an underestimate for the reasons described in Sec. 2.2; such a 
planet would most likely be found around Proxima Centauri or 
Barnard's Star. The probability increases to > 50% when con- 
sidering planets with R P = 1 — 4/?©; these temperate planets 
also most likely to be found around nearby M dwarfs. 

3.2. Ability to Quantify the Local Exoplanetary Demographics 

Several assumptions were made in Sec. 2.2 to extrapolate the 
Kepler population statistics of Howard et al. (2012) to longer or- 
bital periods and cooler stars, and in this section I test the effects 
of these assumptions. The results of these tests are all listed in 
the second section of Table 1 . 

Howard et al. (2012) find roughly equal occurrence frequen- 
cies for 1 — 2/?© and 2 — 4/?© planets. The baseline scenario 
described in Sec. 3.1 conservatively assumes that 2 — 4/?© plan- 
ets and 1 — 2R® planets occur with equal frequency, an as- 
sumption also consistent with more recent results (Dong & Zhu, 
2012; Fressin et al., 2013). For larger radii, Howard et al. (2012) 
find a planet frequency xRp 2 (around sun-like stars). It now 



Perhaps the least conservative assumption in my baseline 
scenario is the decision to assign planets to stars with T e e < 
3600 K, below the range explored by Howard et al. (2012). 
Large-scale radial velocity campaigns have been conducted for 
a few such stars and have found no planets down to ~ 10M© in 
short periods (e.g. Proxima Centauri and Barnard's Star; Endl 
& Kurster, 2008; Zechmeister et al., 2009; Choi et al., 2012). 
Table 1 shows that even if all cool stars have zero planets, 
roughly five planets would still remain for future studies. 

The immediately preceding assumption is of course overly 
pessimistic, since planets with m sin i > 10 M© have been dis- 
covered around a number of M dwarfs (e.g., Butler et al., 2004; 
Bonfils et al., 2005, 2007, 2012). Furthermore, a systematic ra- 
dial velocity survey of a large number of cool (r e jj < 3600 K) 
has recently indicated a planet occurrence frequency of roughly 
90% form sin i= 1-100M©, P < 100 d (Bonfils et al., 2013). 
Using the planet population statistics reported in Table 1 1 of that 
survey (instead of extending the Howard et al., 2012 frequencies) 
for stars with T e g < 3600 K, I predict numbers and types of 
planets amenable to high-contrast characterization nearly identi- 
cal to that predicted by the baseline scenario (see Table 1). In this 
simulation I assume log-uniform frequencies of 2.7% and 3% for 
planets with 10 2 — 10 3 and 10 — 10 2 M©, respectively, while the 
least massive planets (1 — 10 M©) have frequencies of 36% for 
P = 1 — 10 d and 52% for longer periods. I convert m sin i to Rp 
as in Sec. 3.1 (Lissauer et al., 201 1), and impose a maximum ra- 
dius of 1 .3 Rj. The close agreement between this simulation and 
the baseline scenario of Sec. 3.1 hints that the results of Howard 
et al. (2012) and Bonfils et al. (2013) are consistent for small 
planets orbiting low-mass stars under fairly simple assumptions. 
A conclusive test of such a claim is, however, beyond the scope 
of this work. 
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3.3. Effects of Instrumental Capabilities 

In this section I examine how the accessible exoplanet popula- 
tion depends on various observational parameters: telescope di- 
ameter D, inner working angle (IWA), and observing wavelength 
A. The results of all these tests are listed in the bottom section of 
Table 1. 

3.3.1. Angular Resolution: Diameter and Inner Working 
Angle 

As discussed in Sec. 2.3, several GSMT designs have been pro- 
posed, each with different diameters D. The baseline instrument 
used in this work has an IWA defined as 3A/D, so the angu- 
lar magnitude of this IWA decreases as D increases. Planets at 
the smallest separations have the highest observable contrast (in 
both scattered light and thermal emission), so naturally Table 1 
shows that the number of accessible planets increases (by nearly 
a factor of two) when D increases from 30 m (TMT) to 39 m 
(E-ELT). The number decreases by a similar factor when D is 
reduced to 25 m (GSMT). 

Table 1 also shows that at most a single short-period planet 
is predicted to be accessible to high-contrast instruments on cur- 
rent 8 m telescopes, //instruments on these telescopes can reach 
contrast levels of 5 x 10 -9 . In this case, the work presented here 
predicts a > 1% chance of accessible planets for only four stars: 
a Cen AB (15-20% each), Procyon A (5%), and Sirius A (3%). 
These expectation values drop by only a factor of 2 if such fa- 
cilities can only achieve contrasts of 10~ 8 . These star systems 
should therefore be high-priority targets for deep observations 
in imminent 8 m high-contrast surveys. 

Similarly, using a narrower or wider IWA than the baseline 
scenario of Sec. 3.1 respectively increases or decreases the num- 
ber of accessible planets, as expected (see Table 1). However, 
varying IWA has a weaker impact than varying D, even when 
the angular scale of the IWA is held constant. This result may at 
first seem surprising: 2A/D on a 30 m telescope gives a smaller 
angular separation than does 3A/D on a 39 m telescope, yet more 
planets are accessible for characterization under the latter con- 
ditions than under the former. The increased sensitivity at larger 
D for constant angular IWA results from the D 1 1 dependence of 
PSF speckle amplitude on phase and amplitude errors (Guyon, 
2005). Thus for two high-contrast instruments with equal IWA 
(in angular units), the instrument with the larger primary mirror 
can reach better contrast levels. 

3.3.2. Observing Wavelength: Scattered Starlight and 
Thermal Emission 

The choice of observing wavelength A has a number of impor- 
tant consequences: with D it sets the angular scale of the IWA 
(typically specified in A/D), it determines the stellar flux avail- 
able for sensing and correction of phase and amplitude errors, 
and it modifies the planet/star flux contrast (in this work, by de- 
termining the planet's thermal blackbody emission). Naturally, 
observations at multiple infrared wavelengths are also desirable 
to observe spectral features of different molecular species and/or 
surface compositions (e.g., Kuiper et al., 1947; Hu et al., 2012). 

Table 1 shows that the number of planets accessible to high- 
contrast characterization drops steadily as the observing wave- 
length increases from 1 .0 ymi to 3.5 pmi. The detected planets are 
cool and seen in reflective light, so their planet/star contrast re- 
mains approximately constant even as angular resolution and the 
photon flux available for AO wavefront sensing and correction 



both decrease. Nonetheless, Table 1 suggests that several planets 
should exist which are accessible to characterization across from 
at least 1 — 2.2 fim. The number of Earth/Venus analogues ex- 
pected to accessible to characterization is a maximum at 1 .6 yon, 
a value which presumably represents a tradeoff between angu- 
lar resolution and the stellar flux available for wavefront sensing 
and correction. 

Intriguingly, thermal-infrared observations appear to be 
more attractive than observations at intermediate wavelengths 
such as 3.5 fim. All three GSMT projects have proposed AO- 
assisted MIR instruments (Brandl et al., 2010; Tokunaga et al., 
2010; Okamoto et al., 2010; Hinz et al., 2012). The longer 
wavelengths at which these instruments would operate roughly 
correspond to the Wien peak of cool planets' blackbody emis- 
sion, thereby providing a method to measure radiometric radii 
(Morrison, 1973), the surface compositions of airless planets 
(Hu et al., 2012), and the potential to detect exoplanetary ozone 
(e.g., Rugheimer et al., 2012, and references therein). 

The MIR instrument simulated here is not conservative in 
one respect: as in this work's other simulations, I assume that the 
instrument's wavefront sensor and science instrument operate at 
the same wavelength. Current designs for such instruments typi- 
cally rely on facility AO systems, whose wavefront sensors oper- 
ate at shorter wavelengths. Under such a strategy, chromatic ef- 
fects may degrade the achievable contrast to the point that short- 
period planets become inaccessible (Guyon, 2005). Instrument 
designers should therefore consider the feasibility of MIR wave- 
front sensing for these instruments, either with a separate MIR 
wavefront sensor or (to minimize non-common wavefront er- 
rors) via phase retrieval and correction using the main science 
camera (e.g., Malbet et al., 1994; Borde & Traub, 2006). 

Fig. 4 shows the planet population expected around r Cet and 
the expected sensitivity limits for an E-ELT/METIS-like instru- 
ment (A = 10 fim, D = 39 m, and a conservative minimum con- 
trast of 10~ 7 ; Brandl et al., 2012). Note that the planet candidate 
t Cet f (Tuomi et al., 2012) would be amenable to high-contrast 
characterization at 10 fim, assuming a reasonable albedo (< 0.7). 

A high-contrast MIR instrument has a much-reduced angular 
resolution compared to NIR observations, but planetary thermal 
emission boosts the planet/star contrast ratios significantly above 
the values seen at shorter wavelengths; consequently, compara- 
ble numbers of planets are expected to be accessible to NIR and 
MIR characterization. Fig. 5 shows the accessible local planet 
population as a function of Rp and T eq for the thermal imaging 
simulation, and the last row of Table 1 summarizes these results. 

High-contrast thermal-IR instruments thus offer the poten- 
tial to directly measure thermal emission from short -period plan- 
ets down to very small planetary radii. Thermal imaging shifts 
the distribution of sampled planets to substantially cooler T eq 
(< 200 K) as compared to the NIR baseline scenario of Sec. 3.1 
(which is sensitive to rather hotter planets at smaller separa- 
tions). The number of 1 — 2R<$ planets accessible to thermal 
infrared characterization is comparable to the number expected 
in the NIR baseline scenario. However, even these conservative 
planetary population statistics imply a 40% chance of finding an 
Earth/Venus analogue (R P = 1 - 2R®, T eq = 200 - 250 K); 
this probability derives mainly from K-type stars with T s s in the 
range 3900-5500 K, and so may be less sensitive to the assump- 
tions made in Sec. 2.2 for planet frequencies around cooler stars. 
As discussed in Sees. 2.2 and 3.2, if 1 — 2R® planets are 4x 
more frequent than 2 — 4 R® planets then the expectation value 
increases to 1.6. Thus it may be considered likely that such an 
object will be found; such an event would excite considerable 
interest. 
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Fig. 4 Same as Fig. 1 but for t Cet in the thermal imaging sce- 
nario described in Sec. 3.3.2, observing at lOyum with an instru- 
ment similar to E-ELT/METIS. The planets here are detected 
mainly via their thermal emission. From left to right, the star 
symbols and error bars show the contrast expected for planets 
candidates b, c, d, e, and f (Tuomi et al., 2012) over the range 
of albedos and recirculation parameters described in Sec. 2.4. If 
confirmed, planet f would be amenable to direct characterization 
of its thermal emission. 
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Fig. 5 Predicted number of planets (as a function of Rp and r eq ) 
expected to be characterized for the 10 fim imaging scenario de- 
scribed in Sec. 3.3.2. The colors and the small numbers in each 
cell both indicate the expectation value for that combination of 
Rp and T eq , and the histograms show the marginalized distri- 
butions. Compared to the planet population accessible in the 
baseline NIR case shown in Fig. 3, thermal infrared observa- 
tions mainly detect thermal emission from cooler planets around 
warmer stars. Note that the scenario simulated here conserva- 
tively assumes that planets in the two smallest Rp ranges occur 
with equal frequency (see Sees. 2.2 and 3.2). 



4. Conclusions and Future Work 

These simulation of a mostly generic high-contrast instrument 
indicate that, with reasonable extrapolations of planet frequency 
rates from the Kepler mission and/or radial velocity surveys 
(Howard et al., 2012; Bonfils et al., 2013), of order 10 short- 
period planets will be accessible to characterization via high- 
contrast imaging and spectroscopy from the next generation of 
ground-based giant segmented-mirror telescopes (GSMTs). In 
Table 1, the accessible planet population is summarized and 
shown to be insensitive to moderate changes in the assumptions 
made about planet frequencies (see Sec. 3.2) and future instru- 
mental capabilities (see Sec. 3.3). 

Targets for future high-contrast characterization are already 
known today. At wavelengths near 1.2 fim, the planets GJ 876b 
and c (see Fig. 2) and planet candidates a Cen Bb (see Fig. 1) and 
t Cet b should be detectable to the NIR baseline instrument con- 
sidered in Sec. 3.1. GJ 139c and d and r Cet c and d should also 
be detectable if their geometric albedos are > 0.15 — 0.4. Thus 
perhaps a half-dozen suitable planets with a range of masses 
(and presumably radii) are already known; this number is large 
enough that one suspects the modest predictions presented in 
this work may underestimate the true size of the accessible exo- 
planet population. This could be the case if the local exoplanet 
population differs substantially from that observed by Kepler, 
or if the assumptions made in Sec. 2.2 are shown to be invalid. 
Future high-precision radial velocity observations, especially of 
M dwarfs, will enable further tests of these predictions. 

Thermal infrared observations will be desirable, because 
they offer the possibility to directly characterize the thermal 
emission spectra and energy budgets of these short-period plan- 
ets, and thus (via radiometric radii) break the radius-albedo de- 
generacy inherent in reflected light observations. In particular, 
thermal infrared observations significantly improve the chances 
of finding Earth or Venus analogues (planets withRp = 1 —2R® 
and r eq = 200 — 250 K; see Table 1, and cf. Figs. 3 and 5). The 
planet candidate t Cet f (see Fig. 4) is the only currently known 
planet which could be characterized by these mid-infrared obser- 
vations, and radial velocity surveys seem likely to provide addi- 
tional targets in the coming decade. 

The opportunity to observe and characterize the atmospheres 
of relatively small (R P < 4/?©) and cool (T eq < 300 K) plan- 
ets is an exciting one, not least because such planets could be 
similar to Earth and are particularly compelling objects in the 
study of exoplanetary habitability. A useful avenue for future 
study will be to examine in more detail what discernible dif- 
ferences might be used to distinguish Earth-like and Venus-like 
planets (which exhibit similar infrared emission spectra) on the 
basis of mid-infrared observations. One initial avenue might be 
to extend planet classification efforts using optical colors (Crow 
et al., 201 1; Hegde & Kaltenegger, 2012) to longer wavelengths. 

Future, more detailed simulations are required to estimate 
the signal to noise with which the accessible planets discussed 
here can be characterized. Kawahara et al. (2012) have presented 
an excellent first step toward this goal with their recent deter- 
mination that GSMTs will allow low-resolution spectroscopy of 
small planets at 1.27 fim in just a few hours. Further such work at 
both shorter and longer wavelengths will help determine whether 
spectroscopy will be possible, or whether characterization will 
be limited to broadband photometry. Such estimates would also 
permit predictions of the precision with which high-contrast ob- 
servations could measure short-period planets' orbital properties 
(via astrometry), radiometric radii (via MIR observations), and 
perhaps eventually measure rotation periods (Ford et al., 2001), 
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variations in climate (Cowan et al., 2012), and perhaps exoplane- 
tary satellites (Robinson, 201 1; Gomez-Leal et al., 2012). Future 
efforts would also be useful to extend the work presented here 
to observations in polarized light (Zugger et al., 2010). Finally, 
another interesting project would be an assessment of the de- 
tectability of extrasolar comets (Jura, 2005). 

By the end of the current decade, JWST may obtain trans- 
mission spectra of a few temperate, 1-3 /?© planets around M 
dwarfs (Deming et al., 2009; Kaltenegger & Traub, 2009), but 
because transmission spectroscopy probes the highest pressures 
such observations cannot probe deeply into the atmosphere - es- 
pecially if opaque, high-altitude hazes are as ubiquitous as they 
seem to be in the atmospheres of hotter planets (Bean et al., 
2011; Berta et al., 2012; Pont et al., 2012; Benneke & Seager, 
2012). Nonetheless, occultation measurements will be essential 
to inform models of short-period planetary atmospheres and to 
place future high-contrast GSMT observations in a broader con- 
text. 

Future observational efforts are necessary to extend the 
planet demographics measured by Kepler's initial planet can- 
didate sample (Howard et al., 2012) to smaller planetary radii, 
cooler stars, and (less importantly for this paper) longer or- 
bital periods, and indeed such efforts are already being pub- 
lished (Dong & Zhu, 2012; Fressin et al., 2013). Focused high- 
precision radial velocity campaigns of the type used to discover 
planet candidates around a Cen B and t Cet (Dumusque et al., 
2012; Tuomi et al., 2012) should be undertaken for the cool 
stars nearest the sun, both to compare the frequency of plan- 
ets in the solar neighborhood to that observed by Kepler, and to 
discover additional planets for high-contrast GSMT characteri- 
zation. High-precision spectra of all stars within 8 pc of the Sun 
should also be obtained to characterize their stellar properties, 
and thus better estimate the likelihood that they host planets. 

Planned GSMTs will have the potential to measure NIR albe- 
dos, MIR emission, and consequently the molecular chemistry, 
energy budgets, and radii, of significant numbers of small, tem- 
perate, and potentially Earthlike planets around both M dwarfs 
and solar-type stars. Given these positive prospects, continued 
efforts toward a more comprehensive theoretical understanding 
of temperate exoplanet atmospheres is certainly warranted. 
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Crossfield: GSMT High-Contrast Characterization of Short-Period Exoplanets 



Table 1 . Number of Short-Period Planets Accessible to High-Contrast Characterization 



Simulation Parameters 


Planet Size, R P /R B 
1-2 2-4 4-8 


Equilibrium Temperature, r cq /K 
< 200 200 - 300 > 300 


Neva a Total b 


Baseline Scenario (see Sec. 3.1): 

Effects of Astrophysical Parameters (see Sec. 3.2 

Many Small Planets (/i_ MsB = 4/ 2 _ 4 « e ) 

More Longer-Period Planets (for P > 50 d) 

No Planets Around Cool Stars (r eff < 3600 K) 

RV Statistics for Cool Stars (T eff < 3600 K) 


1 + 1 5+ 2 1+ 2 

4 + 2 5+3 it] 
1 + 1 5 + 2 2 + _] 
1+1 3+2 1+1 

i+; a + _[ i±i 


2+2 2+y 4 + 2 

3+, 3+ 2 6+ 2 
4 + 2 3+2 4 + 2 
0+o 1 + 1 4+ 2 
2+2 2+1 4+^ 


0.06 

0.24 
0.06 
zero 
0.05 


8 + 3 

11+t 
10+] 
5 + 2 
8 + 3 



Effects of Instrumental Parameters (see Sec. 3.3): 



Larger Diameter (39 m) 


2+ 2 


9 + 3 


2t] 


4 + 2 


4 + 2 


7+ 2 

'-3 


0.14 


15 + 4 


Smaller Diameter (25 m) 


0+ 1 


3+ 1 


1 ± 1 


1^ 


1 ± 1 


2+ 2 


0.04 


5 + 2 


Much Smaller Diameter (8 m) 


















zero 


0+ 1 


Narrower IWA (2 A/D) 


1^ 




2 +i 

z -2 


3+i 


3+ 1 

—2 


6+ 2 


0.07 


11+ 4 

11 -3 


Wider IWA (4 A/D) 


1+? 


3+ 2 


1 + 1 


2t\ 




3+ 1 

— 2 


0.06 


6_ 2 


Shorter Wavelength (l.O/rni) 


1 ± 1 


5 + 2 




2t\ 




4+3 
H -2 


0.05 


9 + 3 


Longer Wavelength (1.6/jm) 


1 + 1 


4 + 2 


i ± i 


2+ 2 


2+ 1 

^-2 


3+ 1 


0.08 


7 + 3 


Longer Wavelength (2.2 \mi) 


li? 


2+ 2 


i + 1 


2+ 1 

— 2 


1 + 1 


2+ 1 


0.05 


5 + 2 


Longer Wavelength (3.5 /im) 


1+? 


^] 


i« 


l ± l 


1+? 


2+ 2 


zero 


3+ 2 


Thermal IR, Larger Diameter (lOyum, 39 m, 10 -7 ) 




3 + 2 


i ± i 


3 + 2 


1 ± 1 


1 ± 1 


0.41 


6_ 3 



"Expected number of accessible Earth/Venus analogues (R P = 1 — 2R®, T c[| = 200 — 250 K), assuming planets with R P = 1 — 2i?© 
and 2 — 4i?e occur with the same frequency (see Sees. 2.2 and 3.2). 

b Because of rounding, the sum of each row's values may not match the Total value listed. 
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